Non-Markovian behavior of small and large complex quantum systems 
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The channel induced by a complex system interacting strongly with a qubit is calculated exactly 
under the assumption of randomness of its eigenvectors. The resulting channel is represented as an 
isotropic time dependent oscillation of the Bloch ball, leading to non-Markovian behavior, even in 
the limit of infinite environments. Two contributions are identified: one due to the density of states 
and the other due to correlations in the spectrum. Prototype examples, one for chaotic and the 
other for regular dynamics are explored. 
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Introduction- Complex quantum systems are of 
paramount importance in the description of correlated 
many-body systems, such as the ones encountered in con- 
densed matter, as well as few or single body chaotic sys- 
tems. The exact description of such complex systems is 
often not possible because it is cither unfeasible due to 
many degrees of freedom involved, or impossible because 
we do not know all the details of the microscopic model. 
Frequently we are also interested only in the dynamics 
of few degrees of freedom within a larger system. Un- 
fortunately though, even in this case exact solutions are 
very rare. Under certain conditions, which are fulfilled 
in many important situations, one can use approximate 
methods. Such is the case if the central system of in- 
terest is only weakly coupled to the environment with 
fast decaying correlations. This leads to the descrip- 
tion with a relatively simple Markovian Lindblad mas- 
ter equation p], implying a system without memory in 
which information flows only out from the central system. 
While specific models are known in which the reduced dy- 
namics is not Markovian, general understanding is still 
lacking. Such questions resulted in a flurry of recent stud- 
ies of non-Markovian behavior [2]-[5] and characterization 
of reduced dynamics in general [6j [7] . 

In the present work we shall derive an exact description 
of the reduced dynamics of a single qubit immersed in a 
complex system, undergoing unitary evolution. Our goal 
is to characterize the one-qubit channel induced by this 
unitary evolution. We shall assume that the eigenvectors 
of the Hamiltonian governing such evolution can be well 
described by a random unitary matrix. This is a very 
good approximation if the system is quantum chaotic [8] , 
but is also valid under more general circumstances. 

Our main result can be expressed in a very simple ge- 
ometrical picture. The derived one-qubit channel can 
be imagined as an isotropic shrinking of the Bloch ball. 
The radius of this Bloch ball however does not decrease 



monotonically with time but instead oscillates, causing 
non-Markovian behavior. The oscillations are due to (i) 
diffraction on the spectral density and, (ii) due to corre- 
lations between eigenenergy levels. Surprisingly, the first 
contribution will in general lead to non-Markovian behav- 
ior even for an infinite environment. Comparing the con- 
tribution due to eigenenergy correlations leads us to con- 
clude that in the setting studied, chaotic systems display 
stronger non-Markovian behavior than regular ones, as 
quantified by measures proposed in [3J 0] . We also show, 
via exact expressions, that the channel is self-averaging 
for large sizes, meaning that non-Markovian behavior can 
be observed in individual system instances. 

Setting- We study a system of dimension N , divided 
into a single qubit and the rest, acting as an environment 
to which the qubit is strongly coupled. The evolution 
of the total system is determined by a Hamiltonian H. 
The only requirement on H is that the statistical prop- 
erties of its eigenvectors are described by a random uni- 
tary matrix, which is connected to a maximum entropy 
principle |15j . This is conjectured to happen for chaotic 
systems in the semiclassical limit, and is true, by con- 
struction, for the random matrix ensembles [9] suitable 
for describing statistical properties of quantum chaotic 
systems [5]. In quantum information language we want 
to characterize the quantum channel acting on the qubit. 
Once this is done we can study, for instance, whether the 
channel is markovian or not. 

Assume that the initial state of the system is a factor- 
izable state, with a projector in the environment. Other 
choices of initial states will be discussed later. The state 
at later times is thus simply 
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where J7* = cxp (— iHt) (we set h = 1). This induces a 
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where i,j = 0, . . . ,3 with a 3 = {a K ,a y , cr z , 1}. 

Analytic derivation- We are interested in obtaining ex- 
plicit expressions for eq. pj). Writing H in its eigenbasis 
as H = W &id,g{Ei)W\ where W is the unitary matrix of 
eigenvectors of H, we are interested in properties of A(') 
for a unitarily invariant ensemble of Hamiltonians where 
W is a random unitary matrix. We shall calculate the av- 
erage values of all matrix elements of channel A w as well 
as its fluctuations. One finds that given the invariance of 
H under unitary rotations, the average channel, ie., af- 
ter averaging over the unitarily invariant Haar measure 
of W, denoted by must acquire a diagonal form in 
the Pauli basis (which can also be checked by an explicit 
calculation). Such channel is called depolarizing channel 
in quantum information. The matrix (Aj k )u is therefore 
diagonal with time-dependent elements 



a(t) := (Agh u = (4{) u = <A«>„. 



(3) 



Trace preservation means that (A^ 3 )u = 1 and 

(Ag*j =0 j 2 }u — 0. All the physical information about 
the average channel, like the presence of non-Markovian 
effects, is contained in a(t), which is the radius of 
the evolved Bloch sphere of the qubit. The calcula- 
tion of a(t) proceeds by separating the dependence of 
[7* on the spectra Ei and the eigenbasis W, £/* = 
VFdiagexp(— iEit)W\ to obtain 

a(t) = e-^-^iWo^W^W^W^j 

+ Wo», i W? 0ti W^ j W 00 , j ) w 

where we have used both Einstein's summation conven- 
tion and tensorial notation. Latin indices run over the 
whole system, whereas Greek ones run over the system 
minus the qubit. One can then average over the uni- 
tary Haar measure of W using the exact formulas in [10] , 
obtaining the exact expression 



a(t) 



N 2 \f{t){< 
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(4) 



with f(t) — i J^j ex P ( — iEjt) being the Fourier trans- 
form of the level density. The details of the calculation 
are to be found in the additional material |12) . 

The evaluation of the fluctuations of matrix elements 
is of interest, as it indicates how a single member of the 
ensemble will resemble the behavior of the ensemble av- 
erage. Its calcultion involves 8-point correlations of W, 
and the Weingarten function for permutations on 4 ele- 
ments, which we have calculated [TT]. Let us define by 
{[^j\] 2 )u — tne standard deviation of ma- 



trix element A^'jj,. Again, due to the symmetry there are 



only three different fluctuations: those of diagonal ma- 
trix elements, those of off-diagonal elements in a 3 x 3 
block Aj k and those of A-*g. The exact expressions to 
all orders in 1/ A'' is given in the additional material |12j . 
here we only give the leading terms in which are 



^2 2 
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i + \f(t)\ 4 -(f*(trf(2t)+f(tfr(2t)) 
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with i,j = 0,1,2. Equations Q and ^ constitute our 
main result. 

In the above results we have taken the initial state of 
the environment to be a projector. Due to the unitary 
invariance we can choose for IVWv) (VWvl any state. Be- 
cause (Aj k )u is linear in the initial state, any convex 
sum of projectors, ie., a density matrix of the environ- 
ment, will also lead to the same average channel. Fluc- 
tuations though, which are not linear in the initial state, 
do change. In particular, the size of the fluctuations will 
scale as ~ l/(Nr), if r is the rank of the initial state 
of the environment. For instance, if the initial state of 
the environment is an identity matrix, corresponding to 
the environment at high temperature, the fluctuations 
scale as ~ 1/A^ 2 instead of ~ 1/N as for the projector, 
meaning that self-averaging is stronger. 

A random matrix example- We illustrate the above 
results by taking H from the Gaussian Unitary Ensemble 
(GUE). This kind of Hamiltonians have been success- 
fully used to describe a wide range of physical systems 
including chaotic systems, condensed matter systems 
and quantum environments [HI O H3]- For t — ¥ oo the 
induced channel AW is closely related to the so-called 
random quantum channel, in which {/' is replaced by a 
random unitary. Random quantum channels are usefull 
in quantum information theory |16j and have been 
used to prove that the conjecture about superaditivity 
of channel capacities is false [T7]. Because the joint 
probability distribution of eigenvalues is known for GUE 
we can perform explicit averaging over the spectrum, 
obtaining an expression for the average f(t); note that 
due to self-averaging for large N the average behavior is 
observed also in individual samples. As Aw is quadratic 
in f{t) it can be expressed in terms of 1- and 2-point cor- 
relations, which are known exactly for any dimension [5]. 
Strength of the interaction is fixed by (\Hi ,j\ 2 ) = j?, 
resulting in the spectral span of 4 (determining the 
shortest time scale) and the Heisenberg time being 2A^ 
(giving the longest time scale, i.e., the inverse level 
spacing) . The level density is i?i (E) = J^f^Q 1 <Pj C®) > 
where wAx) = ?r Nx2/l %,( x ^jN/2) and H< 

^2ij]y/2n/N 

are Hermite polynomials. The cluster function, 
giving correlations between different levels, is 
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FIG. 1: (a) Theoretical dependence of a(i) for Poisson ex- 
ample, Eq.Q for N = 4, 8 and oo. (b) Same three sizes for 
GUE example, Eq.|6|. We also show three diagonal elements 
of A (t) for one instance of N = 12000 (three thin blue curves) 
and theoretial fluctuations, displayed around N — oo theory 
as a gray shadow. 



for GUE T 2 (E l ,E 2 ) = (j^f^ ^(^1)^(^2)) • 

One can show that iV 2 (| /(i)| 2 ) G UE = N + 
J dE 1 dE 2 e- i ^- E ^ t [R 1 (E 1 )R 1 (E 2 ) - T 2 (E U E 2 )}, 
which can be evaluated explicitly for any N. Let us 
define b 2 (t) = (l/N) J dE 1 dE 2 e- i< - El - E ^ t T 2 (E 1 , E 2 ) 
and foi(i) = (l/N) J dEe~ iEt R^E) . Normalization is 
such that 6i(0) = 1 and 62(00) = 0. The final formula is 



<«(*)> 



GUE 



N 2 b\(t)+N(l-b 2 (t))-l 
N 2 -l 



(6) 



Each of the contributions approach a simple expression 
in the limit N — > 00: limjv^co bi(t) = Jl ^ , while the 
leading order of the form factor 1 — b 2 (t) is t/2N for 
t < 2N and 1 otherwise. 

We can see that there are two contributions to 
(a(i))GUE- The first one comes from the Fourier trans- 
formation of the energy density. The second one, given 
by the form factor and being due to cigcnencrgy corre- 
lations, is of the order l/N compared to the first one, 
and is therefore important for moderate TV. For large 
N the second term can be neglected in Eq. [6j giving 
lim A r^ 00 (a(t)) GUE = [Ji(2i)/t] 2 . The form of (a(i))cuE 
for small and large N in shown in Fig. [TJd . There we also 
show fluctuations, which can be for large N obtained by 
using f(t) ps Ji(2t)/t in Eq. [5] . The fluctuations de- 
cay with the system size as ~ l/N. Therefore, for suf- 
ficiently large system the fluctuations are smaller than 
( Q! (0)gue, i-e., the dynamics is self-averaging. Even tak- 
ing a single member H of the GUE ensemble one gets 
the average behavior (a(i))GUE, as can be seen in Fig. [I] 



for N — 12000. If we would take the initial state of the 
environment to be the maximally mixed state, instead 
of |'f/ , env)('i/'env|j similar self-averaging would be achieved 
already for N s» -\/l2000 rj 100, which is about 7 qubits. 

Poisson example- As a second example we show one 
still possessing unitary invariance, but having Poissonian 
eigenenergies with no correlations [2], and with a flat 
level density, being a model for regular systems. The cal- 
culation goes exactly as in the previous example. Taking 
into account that there are no correlations among differ- 
ent levels and the spectral density is flat (T 2 (Ei, E 2 ) = 0, 
R X (E) = NQ (\E - 2|) /4, where 9 is the Heaviside step 
function), we get (see Fig.[lji) 
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Non-Markovian behavior- Having calculated a(t), 
one can immediately draw conclusions about the non- 
Markovian behavior of the channels. Consider the map 
that takes a state from a time t to t + r, A( 4, * +T ) = 
A (*+-) [A«] . This is, in general, not a physical map, 
which implies that the trace one operator associated via 
the Jamiolkowski isomorphism J is not a physical state. 
In [3], the deviation of positivity for such operators is 
taken as a measure of non-Markovian behavior Ai 1 . We 
define 



](t) = lim 

£->0+ 



\\J(A {t+e 



- 1 




if a(t) > 
otherwise 

(8) 

which will be positive whenever a(t) increases (the de- 
tails are presented in the supplementary material [12]). 
With this figure of merit one can calculate the values 
of Ali = J °° g(t)dt. A different criterion is based on 
the evolution of distinguishability of states with time [5] 



and is defined as M2 = max 



Xr>o dicr ^°,Pi,i), 



P0,l(°) Ja>0 

where a(po(t), pi(t)) is derivative of the trace distance 
between po,i(t). The states po,i that maximize such 
quantity for our channel are any two orthogonal pure 
say pi — \i)(i\. In such case pi(t) = ^ ±a ^ a * 
The last measure to be ex- 



states, 



and M 2 = 2 / 4>0 dtd(t). 
amined quantifies non-Markovian behavior via the non- 
monotonicity of entanglement decay of our qubit with 
an ancilla qubit £Q and is as such, as we will see, 
weaker than M\ i2 : M3 = Jc >0 dt[dC(p(t))/dt], where 
C is a measure of entanglement (to be taken here as 
the concurrence [H]), p(0) is a Bell state in the two 
qubits and the quantum channel acts on a single qubit. 
The concurrence for the corresponding state will be in 
our case C(p(a)) = max{0, (3a — l)/2}. The final re- 
sult is M.3 = § / d>0 01/3 <idt. In table U\ we report 
several values of all three measures for different envi- 
ronments. We can see that both A4\ v2 indicate non- 
Markovian behavior exactly at times when a(t) increases, 
in other words, when the Bloch ball expands. If we 



4 



explicitly write Mi — 3/2Y^ In (a(i/)) — In (a(ti)) and 
Mi = 2 a(tf) — a(ti), where both summations are 
over all intervals [ti,tf\ on which a(t) increases, it is also 
easy to understand why the behavior of M\% is different 
with N . Because of the divergence of logarithm at 0, the 
behavior of Mi is dominated by values of a(ij) which 
decrease with N, eventually becoming for N — > oo, 
causing the increase of Mi with N. On the other hand 
Mi is dominated by terms a(tf) that decrease with TV, 
see Fig. [TJd . For Poisson example Mi increases due to a 
trivial Nj (N + 1) prefactor. Looking back at our results 
and the two examples of a GUE and Poissonian ensemble, 
we can see that for small times non-Markovian behavior 
is due to diffraction on the spectral density. Provided the 
spectral span S is finite, there will always be oscillations 
in a(t) on the time-scale causing non-Markovian 

behavior. How fast these oscillations decay with time 
depends on the singularity at the spectral edge - sharper 
features lead to slower decay of oscillations with time. In 
condensed matter systems singularities at spectral edges 
(van Hove singularities) are quite common. Surprisingly, 
non-Markovian behavior is present even for an infinite en- 
vironment, where one would perhaps expect that there is 
no "back-flow of information" from the environment to 
the qubit. For smaller systems the term with the 2-point 
correlations also leads to non-Markovian effects. Indeed, 
for chaotic systems 1 — hi (t) increases with time, lead- 
ing to an additional increase of a(t). This contribution 
occurs on the time-scale of the inverse level spacing. In- 
teresting to note is, that comparing the GUE case, mim- 
icking chaotic systems, with the Poissonian for regular 
dynamics, shown in Fig. [T] , one can conclude that non- 
Markovian effects are stronger in chaotic systems than in 
regular ones. This is yet-another example of a counter 
intuitive behavior of quantum chaotic system. Another 
is their stability, where quantum chaotic systems can be 
less sensitive to perturbations than regular ones [19] ■ 
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N = 4 


4.375 


0.378 





0.555 


0.156 





N = 8 


6.102 


0.236 





1.064 


0.173 





N —> 00 


00 


0.051 





00 


0.195 






TABLE I: Different values of non-Markovian behavior for sev- 
eral environments. Notice how the two measures Mi,2 have 
different tendency for the GUE case, and how M3 can not be 
used to detect non-Markovian behavior in our systems. 



Conclusion- We analytically calculate a quantum chan- 
nel describing the reduced dynamics of a single qubit 
within a larger system. Unitary evolution by unitarily 
invariant Hamiltonian leads to simple diagonal channel 
that can be visualized as an isotropically oscillating Bloch 
ball. The average value of the diagonal matrix element 



has two contributions: (i) one from the Fourier transfor- 
mation of the energy density, and (ii) from correlations 
between eigenenergies. Provided there is some eigenen- 
ergy repulsion, as is the case in quantum chaotic sys- 
tems, the second contribution will lead to semiclassically 
small non-Markovian behavior. This effect is stronger 
for more chaotic systems. The contribution due to en- 
ergy density in general leads to non-Markovian effects 
even in the limit of an infinite environment. We also 
calculate channel fluctuations, showing that the dynam- 
ics is self-averaging for large systems. This means that 
non-Markovian effects should be observable already in 
small individual systems, making it an exciting exper- 
imental challenge. Acknowledgments- Support by the 
Program Pl-0044, the Grant JI-2208 of the Slovenian 
Research Agency, and projects CONACyT 57334 and 
UNAM-PAPIIT INI 17310 are acknowledged. 
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SUPPLEMENTARY MATERIAL 
Fluctuations of A^ 

After straightforward but tedious calculation, we ob- 
tained exact results for all three different fluctuations. 



They can all be expressed in terms of the Fourier trans- 
formation of the level density denoted by f(t) and its 
powers, 



/W-^E^PH^)- ( 9 ) 

j 

Exact expressions (to all orders in 1/N) for fluctuations 
of matrix elements of are the following (in these ex- 
pressions no averaging over eigenenergies is performed 
yet; therefore, any spectrum can be used): 



2iV(l * )(1 9 ) 



(1 - ±) + (2 - £ - ±) jiVl/WI 4 + >(2*)| 2 - ||/W| 2 



+(i--){[rw] 2 /(2t) + [/(<)] 2 r(2<)} 



(10) 



To get the fluctuation, one has to subtract from the pre- unitary and eigenenergies, has to be performed before 
vious expression (A^) 2 (here the averaging over both, squaring). Other two fluctuations are 
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2N(l-&)(l-&) 



12 



(! - Jp) ~ 3 I/WI 4 - ^2 l/( 2i )| 2 + + {ir(t)] 2 f(2t) + [f(t)] 2 f*(2t)} 



N 2 ' 
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(1-Ja) + I/Wl 4 + ^I/( 2 0I 2 -^I/W| 2 - 



a- ^){\r(t)] 2 f{2t) + [f(t)] 2 rm} 



(ii) 

(12) 



Notice that one gets correlations up to 4th order, includ- 
ing those that mix f(t) at different times, making the ex- 
act averaging over eigenenergies, for instance in terms of 
Hermite polynomials for GUE, very cumbersome. In the 
leading order with respect to 1/N one can, often, forget 
about correlations, and interchange averages of powers 
with the powers of averages. Thus, if h t :— (f{t)) spcctrum , 
we obtain for the leading order 



- 2 (A^) = 

-XI) = 
- 2 (A^) = 



1 + (h* t 2 h 2t + h 2 h* t ) - 3\htl 
2N 

l + |/i t | 4 -(^ 2 fe 2t + /i 2 fe* t ) 
2N 

- 2 (A^). 



(13) 



This approximation is indeed valid in both examples ex- 
amined in the main text for large dimensions. For GUE 
ensemble h t = J\(2t)/t, giving very simple expression 
for correlations. They are shown in Fig. [2] . In Fig. [3] 
we also show the values of 6 off-diagonal matrix elements 
Aj*fci j 7^ k = 0, 1, 2, for one GUE instance of dimension 
N = 12000 (the same data as shown in the main text). 
Because the average values of these off-diagonal elements 
is zero (thick line at in the figure), they simply fluc- 
tuate around with the amplitude given by theoretical 
cr 0) i, Eq.( |13[ ), and shown as a gray shadow in the figure. 

For smaller size N — 4, and again GUE ensemble, the- 
oretical expressions for fluctuations ( 10|ll|12 1 are shown 
in Fig. [I]. One can see that the time dependence is quite 
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FIG. 2: Scaled fluctuations of matrix elements of A*-*' for a 
GUE ensemble and large sizes, where one can use ht for f(t). 




FIG. 4: Theoretical formulas for fluctuations given in 
Eqs.| |10|ll|12) for size N = 4. All is for a GUE ensemble. 
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FIG. 3: Values of off-diagonal elements of A^, with j 7^ k = 



0,1,2, for one GUE instance of size N = 12000 (thin blue 
curves). Fluctuation 00,1 is shown as a gray shadow around 
the average at (thick red line). 



complicated. General feature of fluctuations is that they 
are very small for short times, and reach their maximal 
value before the first revival in a(t) (this comes about 
due to the presence of f(2i) term in fluctuations). 



where averaging is over GUE spectrum. Such average 
can be written as 



(s2 e -i&-B,)t\ =N+ 

I GUE 

J &E l &E 2 e- l{E ^ E ^ t (Y, 8(Ei - E l )8{E 2 - E 



GUE 



The quantity to be averaged is the two point correlation 
function: 

R 2 {E l ,E 2 ) = r£ d 5{E l -E i )5{E 2 -E j )\ , (15) 

Wj / QUE 

which can be expressed in terms of the level den- 
sity of states R\{E) and the two level cluster function 
T 2 (Ei, E 2 ): 

R 2 (E 1 ,E 2 ) = R 1 {E 1 )R 1 (E 2 ) - T 2 {E X ,E 2 ), (16) 

for which explicit expressions in terms of Hermite poly- 
nomials exists. In particular, 



JV-l 

Rx{E)= $>f(£), 
3=0 



(17) 



and 



Details for the GUE calculation 



We want to evaluate the quantity 



' JV-l 



T 2 (E l7 E 2 ) = \ Vj{Ei) Vj {E 2 ) I . (18) 
j=o 



with 
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(14) 



-Nx 2 /4 



Hjixy/N/2), (19) 
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and %j being the Hermite polynomials. Explicit expres- 
sions for moderate iVs can be obtained by straightforward 
calculation with the aid of symbolic computational pro- 
gram. The final expression for the desired quantity is 
thus 



(l/(*)| 2 > 
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N + J AE l dE 2 e- i ^- E ^ t R 2 {E l ,E 2 ) 
N~ 2 



(20) 

where R 2 , R\ and T 2 are given in terms of the Hermite 
polynomials. 

In the large N limit, simple expressions are also avail- 
able. The limit of the level density is known as the semi- 
circle law, and yields an ellipse with semi axis determined 
by normalization. Its Fourier transform is a Bessel func- 
tion J\\ 



d£iM£)e- iE < N^t. (21) 



For the second term, the integral in the large N limit 
yields the well known two level form factor for the GUE: 



J dE 1 dE 2 e- i ^- E ^ t T 2 (E 1 ,E 2 ) 
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JV- |t|/2 Ht<2N 

otherwise ' ^ ' 



Measures of non-Markovian behavior 



The depolarizing channel maps a state p — ^+'" cr — > 
T^a {p) = ^~ +a ^ r ' (y . This is precisely the map correspond- 



ing to {A)u- We shall now work in the Choi basis, which 
is, for a single qubit, {|0)(0|, |0)(1|, |1)(0|, |1)(1|}. The 
matrix representation of such a channel, in the afore- 
mentioned basis is 



/ i+q 
' 2 



\-oc \ 
2 » 

a 







a 

o o i±^y 



(23) 



One can also think of the map from a time t\ to t 2 , which 
is not necessarily physical. The matrix representation of 
such a map is 



V 



*2,tl 



\ 



a r 
a r 
\J - a r \ + a r J 



(24) 



with a r = a(t 2 )/a(ti). The associated state, via the 
Jamiolkowski isomorphism J is 



1 



(\+a r 



t 2 ,ti 







(25) 



a r \ 

\-ol t 

a r | + a r J 

with eigenvalues (1 — a r )/A (three times) and (l + 3a r )/4. 
With this one directly arrives to the result. 
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